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We report on two recent multiloop results in QED: (i) the four- loop corrections to the conversion relations 
between the QED charge renormalized in the on-shell and MS schemes; (ii) analytical evaluation of a class of 
asymptotic contributions to the muon anomaly at five-loops. 
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1. Introduction 

The study of the anomalous magnetic moment 
of the muon a M is a long-standing challenge for 
both theory and experiment. It has been mea- 
sured with impressive accuracy at the level of 0.5 
parts per million [I]: a^ p = 116592080(63)-Kr n 
From theory side the anomalous magnetic mo- 
ment has been studied in great detail through the 
computation of higher order corrections (see, e.g. 
reviews [21314] ). These higher order corrections 
to a*' leo are basically classified into three classes: 
pure QED, electroweak and hadronic contribu- 
tion. 

A discussion of the electroweak and hadronic 
corrections can be found in Ref. [5] and references 
therein. Within this work we consider higher or- 
der corrections to the pure QED part. Starting 
from two-loop, diagrams with internal fermion- 
loops can arise, where the fermion-type of the in- 
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tcrnal loop can in general be different from the 
external muon. The one-loop [6] and two-loop 
|7I8I9I10I11I12I13J contributions have been com- 
puted more than 50 years ago. The three-loop 
order has been computed numerically and ana- 
lytically [141151161171181131201 

The complete calculation of hundreds of dia- 
grams contributing to the muon anomaly at four 
loop order is only possible by numerical integra- 
tion, which was performed in a remarkable long- 
term effort by Kinoshita and his collaborators 
[21 22J. Even some numerically important five- 
loop diagrams have been computed by now [23] , 

Accurate numerical calculation of highly di- 
vergent multiloop Fcynman amplitudes is highly 
nontrivial task. The analytical calculation (of 
even a very particular class of diagrams) should 
be quite useful as an independent check in both 
directions! (see e.g. instructive examples in 
Refs. [24125126127121128129] h 

In particular starting from two-loop there 
arise logarithmic contributions of the type 
log(M fl /M e ) ) where M M is the mass of the muon 
and M e the mass of the electron, respectively. In 
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view of the large mass ratio M„/M e ~ 200 one 
can expect, that these logarithms play a domi- 
nant role. The logarithmically enhanced contri- 
butions arise from the insertion of the electron 
vacuum polarization (eVP) into the first order 
muon vertex diagram, but they can also appear 
through light-by-light (LBL) scattering diagrams. 
Examples for both diagram types are shown in 

Fig. m 



Here 





Figure 1. Example diagrams leading to domi- 
nant logarithmic contributions from electron vac- 
uum polarization (VP) insertions and light-by- 
light (LBL) scattering diagrams. 



In what follows we will consider only eVP con- 
tributions to a M which are produced from dia- 
grams like Fig 1 (eVP) but with the photon prop- 
agator receiving all possible QED perturbative 
corrections made of electron loops and photon 
exchanges only. The general structure of eVP 
was throughly studied long ago in the pioneer- 
ing publications of B. Lautrup and E. De Rafael 
and R. Barbieri and E. Remiddi [30|31j . 

There it was found, that the asymptotic part of 
the muon anomaly a°f , which contains these log- 
arithmic contributions originating form the elec- 
tron vacuum polarization function insertions and 
the mass independent term, can be obtained with 
the help of the electron vacuum polarization func- 
tion in the asymptotic limit M e — » 0. The corre- 
sponding master formula reads 
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In the above formulas II is the photon polariza- 
tion operator (in QED incorporating exactly one 
fermion field), M is the on-shell fermion mass, 
and the fine structure constant a is defined in 
the classical OS-scheme by the condition: 



Il(-q 2 =0,M 2 ,a) = 0. 



(4) 



Once II as is computed to some order in a the 
master formula (fTJ) immediately delivers the cor- 
responding eVP contribution to a M . 

At present this technique has been applied in 
order to find the complete eVP contributions to 
a M up to three-loop and four-loop order [301311251 

m\. 

In our talk we extend these results by one order. 
As a by-product we will also derive the four-loop 
corrections to the conversion relations between 
the QED charge renormalized in the on-shell and 
MS schemes. 

2. Photon polarization operator in MS- 

scheme 

Thus, the main problem is to compute the 
asymptotic photon polarization operator {Q 2 = 



n as (Q/M,a) =5^IP"'W(Q/M) (- 



(5) 



with every Il as '^'(Q/M) being, in faclQ a polyno- 

mial of order not higher than i in £qm = log(-pr). 
It is very useful to consider first the photon 
polarization operator in MS-scheme written as 
II(— q 2 , to 2 , /i, a) where a = a ([i) and to = 
m MS (fi) are the running coupling constant and 
the fermion mass in the MS-scheme, while \x 
stands for the MS renormalization point. Indeed, 



3 This is a consequence of the Weinberg theorem 32 . 
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the MS polarization function II has a smooth 
massless limit: 



II(Q 2 ,m 2 = 0,a) = ^irV/Q 



?0) 



.(6) 



A use of the fundamental concept of the invariant 
charge [33 34] directly leads to the connection 
between II and II os 



a 



OS 



-7, (7) 



l + n os (Q 2 ,M 2 ,a os ) l + II(Q 2 ,m 2 ,a) 

Eq. allows one (see e.g. [HI]) to relate a MS 
and a os through the MS polarization operator at 
zero momentum transfer 



«/ 2/ 2n {<*(v) 



n(Q 2 = 0,m 2 ,a) = ^n W ( / i 2 /m 2 ) 



.(8) 



Thus, to compute the II as at four loops we 
need to know the MS polarization function II in 
the massless and momentum-less limits as well 
as the relation between on-shell and MS masses 
of a fermion in QED at three loops. The latter, 
fortunately, is available from [35136) . 

2.1. Massless limit of IT 

We have computed II (a 2 /Q 2 ) using the par- 
allel version of FORM [37138139] . The contribut- 
ing diagrams were first generated with the pack- 
age QGRAF [40]. The reduction to master inte- 
grals was performed with the help of an auxil- 
iary integral representation [4T] . The result reads 

(^Q = log(^)): 
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2.2. Momentum-less limit of II 

The limit q 2 — > leads to the evaluation of 
massive tadpole diagrams. Their computation 
has been performed with FORM [37142143] based 
programs. All appearing tadpole diagrams have 



been reduced to master integrals with the help of 
Laporta's algorithm [20 44] . The arising polyno- 
mials in the space-time dimension d = 4— 2 e have 
been simplified with the program FERMAT [45 . 
The remaining master integrals are known ana- 
lytically to sufficient high order in s and have 
been taken from Refs. [46I47I48I49I50I51I52I531 



[54] , The following result for II 
found 
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at Q — was 
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Here l^ m = log(^), C„ — ((n) is Riemann's zeta 
function and a n = Li„(l/2) = J2Zi V(2 l «"). 

3. Conversion formulas for a 

Let us define the conversion factor C„ a , which 
converts the fine structure constant a in the MS- 
scheme into a in OS-scheme: a — C^a &'■ 



^ryrv J- 



^m 



(id 



A use of eq. ([10]) directly leads to our result for 
Cyi (for brevity we skip the lower order expres- 
sions for C — with i = l,2,3 which can be found 
in [55]) 
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Note that in the process of deriving (fT2|) one needs 

also n (l V7rn 2 ) and H ( % 2 /Q 2 ) for * = 1,2,3. 
We have not written the corresponding results ex- 
plicitly to save space; the reader could find them 
e.g. inRef. [27]. 



5. eVP contribution to the muon anomaly 
at fifth order 



Combining eqs. (|1I13[) we arrive to our main 
result: the complete eVP contribution to the 
muon anomaly at order a 5 (below we write down 
only the new five-loop result; we have confirmed 
all previously known eVP contributions to a M at 
three and four loops) 



i>2 
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4. Asymptotic photon polarization opera- 
tor in on-shell scheme 



The combined use of eqs. (|9I11I12[) immedi- 
ately leads us to the the following expression for 
the four-loop contribution to the asymptotic pho- 
ton polarizatiorQ 
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Note that in the process of deriving (flB"]) one needs 

also TT W (^ 2 /m 2 ) and n ( %7Q 2 ) for * = 1,2,3. 
We have not write the corresponding results ex- 
plicitly to save space; the reader could find them 
e.g. in Ref. [27]. 



4 All terms proportional to i in \Yi\ could be derived in 
easy way from renormalization group arguments as ex- 
plained in 1 301 . In addition, the process of constructing <$£ 
one should invert the power series for (l + IT as ), which also 
produces remarkably easy-to- compute factorizable fourth 
order contributions to d^ like ct 4 n a3 '( 2 ) Yl as '^ and so on 
(see e.g. |56| and references therein). However, to include 
power suppressed terms of order (M/Q) n to n os is much 
less trivial even for the factorizable contributions |57I58I . 
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Me - log (M M /M e ). Finally, numerically, 
the result reads: 



as, (5) 



63.481 NS - 1.21429 S i. (16) 



Eq. lfl6|) also displays the decomposition of the 
full result into two pieces, corresponding to non- 
singlet and singlet contributions to the photon po- 
larization operator (see, e.g. [59]). 

6. Conclusion and Acknowledgment 



We hope that our main result (|16p for the com- 
plete fifth order eVP contribution to the muon 
anomaly could be of some use for testing the com- 
plicated numerical simulations like it happened a 
few times in the past. 
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Note added. Very recently, two months after 
the conference, the result for the five-loop eVP 
contribution to the muon anomaly from singlet 
diagrams have been published by T. Aoyama, M. 
Hayakawa, T. Kinoshita, M. Nio and N. Watan- 
abe [50] ■ Their result reads: 

A2 (m^/me) (Set I(j)( e ,e) '■ combined) = 
-1.24726(12). (17) 

One observes a good agreement for the singlet 
case between the numerical result (fTTf and our 
prediction ((T6|) . The remaining 2.6 % difference 
presumably comes from power suppressed correc- 
tions to the asymptotic result. 
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